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We report nuclear quadrupole resonance (NQR) studies on the chain Cu sites of PrBa2Cu40g, a 
quasi-one-dimensional conductor with a nearly quarter-filled band. The nuclear spin-lattice relax- 
ation rate 1 /Ti shows a pronounced peak near 100 K caused by fluctuations of electric field gradient 
(EFG). Similar peak was observed for the spin-echo decay rate I/T2, however, at a different temper- 
ature near 50 K. These results and broadening of the NQR spectrum at low temperatures indicate 
that slow charge fluctuations of either electronic or ionic origin freeze gradually at low temperatures. 



There has been increasing interest in quasi-one- 
dimensional (Q1D) correlated electrons. Theoretical 
studies on generalized Hubbard or t-J models for chains 
and ladders have revealed rich phase diagram associ- 
ated with various instabilities towards Mott localization, 
spin density wave, charge order, and superconductivity 
jy, [§, [§. A variety of phases has been indeed observed 
in organic Bechgaard salts Q] and cuprate ladder ma- 
terials H §. Pr-based cuprates, PrBa 2 Cu 3 C>7 (Prl23) 
and PrBa 2 Cu40s (Prl24), also provide good model sys- 
tems of strongly correlated Q1D electrons. Each of these 
has identical structure to the respective Y-based high 
temperature superconductor, YBa 2 Cu3C>7 (Y123) and 
YBa 2 Cu 4 8 (Y124). The Cu0 2 planes in Pr-based com- 
pounds, however, show antiferromagntic order j?], || and 
their contribution to the optical conductivity revealed 
features of a charge-transfer insulator || [h]]. The in- 
sulating nature of the Cu0 2 planes was explained by 
Fehrenbacher and Rice jy] based on a model of the lo- 
calized hole states made of Pr-4/ and 0-2p„- hybridized 
orbitals. Thus the active element for low energy spin 
and charge excitations in Pr-based cuprates are the Cu-0 
chains, single chains in Prl23 and double (zigzag) chains 
in Prl24. 

These two compounds show contrasting transport be- 
havior. Prl24 is a highly anisotropic Q1D metal showing 
extremely large in-plane resistivity anisotropy p a /pb ~ 
1000 at 4 K |J, H, |J|, where b is along the chain direc- 
tion and a is perpendicular to it within the Cu0 2 plane. 
Although p shows monotonic metallic temperature de- 
pendence, p a exhibits a broad peak near 130 K, indicating 
incoherent transport perpendicular to the chains at high 
temperatures. The angle resolved photoemission spec- 
tra (ARPES) also revealed a clear one-dimensional 
dispersion which crosses the Fermi level near the mo- 
mentum k ~ 7r/4, pointing to a nearly quarter-filled 
conduction band. In Prl23, the spectral weight of the 
chain contribution to the optical conductivity H as well 
as the ARPES spectra jl6j are also consistent with an 
approximately quarter-filled conduction band. However, 
semiconducting behavior of the dc-resistivity, steep sup- 



pression of the optical conductivity below 0.1 eV and van- 
ishing ARPES intensity near the Fermi level all indicate 
the existence of a charge gap. Nuclear magnetic reso- 
nance (NMR) and nuclear quadrupole resonance (NQR) 
experiments in Prl23 ]l7| revealed line broadening and 
anomalies in relaxation rates, which were ascribed to a 
charge ordering instability. 

Apparent absence of static charge order in Prl24 may 
be due to slight deviation from the quarter-filling [H| or 
to the geometrical frustration of intersite Coulomb inter- 
actions in zigzag chains Jlq] , However, dynamic signature 
towards charge order instability is indicated by the opti- 
cal conductivity in Prl24, which splits into a zero-energy 
Drude part with a small weight and a finite energy mode 
centered at 40 meV [^0| . In this letter, we report re- 
sults of NQR experiments for the chain Cu sites in Prl24. 
We found a pronounced peak in the temperature depen- 
dences of the nuclear spin-lattice relaxation rate (1/Ti) 
and spin-echo decay rate (1/T 2 ), as well as broadening 
of the NQR spectra, indicating gradual freezing of fluc- 
tuations of electric field gradient (EFG) with decreasing 
temperature. 

The powder sample of Prl24 was synthesized by solid 
state reaction under high pressure as described in [M . 
Standard spin-echo pulse technique was used to obtain 
NMR spectra. The spin-echo decay rate (1/T 2 ) was ob- 
tained from the spin-echo intensity measured against the 
time separation r between tt/2 and ir pulses. The nu- 
clear spin-lattice relaxation rate (1/Ti) was measured by 
inversion recovery method. 

The NQR spectra for the chain Cu sites in Prl24 are 
presented in Fig. El We observed well resolved resonance 
lines for both 63 Cu and 65 Cu isotopes with the full- width 
at half-maximum (FWHM) of 560 kHz for 63 Cu at 300 K. 
The resonance signal from the planar Cu sites was ob- 
served near 70 MHz at 1.5 K as reported elsewhere [||, 
indicating an antiferromagnetic order of the planer Cu 
spins. The NQR spectrum for the chain Cu sites shows 
little change at the Neel temperature of the planar Cu 
spins (~ 220 K). 

The temperature dependence of 1 /Ti for the chain Cu 
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FIG. 2: The temperature dependence of 1/Tj for 63 Cu and 
65 Cu. The data for the chain 63 Cu sites in Y124 @axe shown 



for comparison. 



FIG. 1: NQR spectra for the chain Cu sites in Prl24. The 
peak near 20.2 MHz (18.7 MHz) is due to 63 Cu ( 65 Cu) nuclei. 
The inset shows the temperature dependence of the line width 
(FWHM) for 63 Cu. 



sites is shown in Fig. || Here 1/T± is defined as the in- 
verse time constant for the exponential recovery of the 
NQR intensity after the inversion pulse, which is three 
times larger than the standard definition if the relax- 
ation is due to magnetic processes. The fitting of the 
recovery curve to an exponential function was satisfac- 
tory over two decades except below 100 K, where the 
distribution of 1 jT\ results in a slightly non-exponential 
recovery curve. Above 200 K, 1/T± is nearly proportional 
to T, similar to the results in Y124 Q. Below 200 K, 
however, \ jT\ deviates significantly from the T-linear be- 
havior and shows a pronounced peak near 100 K. 

A similar peak was observed also for I/T2 shown in 
Fig. H The spin-echo decay curves can be fit well to 
an exponential function exp(— 27-/T2) below 200 K. At 
higher temperatures, the spin-echo decay curve contains 
also a small Gaussian component and the fitting becomes 
less satisfactory. Although I/T2 is almost independent 
of temperature above 100 K, there is a strong peak near 
50 K. Note that this temperature is lower than the peak 
of l/Ti by 50 K. The peak value 1/T 2 ~ 120 ms" 1 is more 
than 5 times larger than the value reported for Y124 

The isotopic ratio of these relaxation rates can be used 
to identify the relaxation process. Both \ jT\ and I/T2 
are larger for 65 Cu than for 63 Cu when the tempera- 
ture is sufficiently higher or lower than the peak (Figs. ^ 
and |3|). Their ratio is close to the ratio of squared nu- 
clear gyromagnetic ratio ( 65 ^/ 63 ^) 2 = 1.148, indicat- 



ing that relaxation is caused by fluctuations of the lo- 
cal magnetic field acting on nuclei. In the temperature 
range near the peak, on the other hand, both relax- 
ation rates are larger for 63 Cu, with the isotopic ratio 
close to the ratio of squared nuclear electric quadrupole 
moments ( 65 Q / 63 Q) 2 — 0.856. Therefore, we conclude 
that the peaks in the relaxation rates are caused by 
fluctuations of the electric field gradient (EFG) tensor, 
V a j3 — d 2 V/dx a dxfj, where V is the electrostatic poten- 
tial at the nuclear position. 

The relaxation rates can then be related to the time 
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FIG. 3: The temperature dependence of I/T2 for 63 Cu and 
65 Cu nuclei. 
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FIG. 4: The temperature dependence of the correlation time 
for EFG fluctuations. The line shows a fit yielding the acti- 
vation energy of 430 K. 

correlation functions of EFG components. A simple in- 
terpretation of temperature dependences of 1/Ti and 
I/T2, in particular the difference in their peak temper- 
atures, is given by a classical model of motional nar- 
rowing |2^, [23|. Let us assume a simple exponential 
form for the correlation functions, (V a p(t)V a f3(0)) — 
(V^g) exp(— t/r c ), where r c is the correlation time of the 
fluctuations. The spin-lattice relaxation rate is given 
by their Fourier transform at the NQR frequency uj n 
(20 MHz in our case), 



where A is the magnitude of the transition matrix el- 
ement of the quadrupolar interaction, which is propor- 
tional to the amplitude of the EFG fluctuations (V a p) 
appropriately averaged over different tensor components. 
If A is constant but t c depends on temperature, 1/Xi at- 
tains the maximum value when uj n T c = 1. Thus the peak 
of 1 /T\ is naturally understood due to increase of r c with 
decreasing temperature, i.e., slowing down of the EFG 
fluctuations at low temperatures. By applying Eq. ([[]) to 
the experimental data, we deduced A/2ir = 170 kHz and 
1/t c is obtained as a function of temperature near 100 K 
as shown in Fig. ^. Although the temperature range is 
limited (75 K > T > 140 K), we can extract an activa- 
tion energy 430 K from the temperature dependence of 
1/t c . 

In contrast, it is known that I/T2 determined from 
the spin-echo decay curve takes the maximum value 
1/T 2 ~ A when Ar c ~ 1 @ §4|. Since A is two 
orders of magnitude smaller than u) n , gradual slowing 
down of the EFG fluctuations accounts for the fact the 
the peak in I/T2 occurs at a lower temperature than 



the peak of 1/T\. The plot in Fig. || indeed shows that 
1/t c ~ 10 6 s _1 ~ A near 50K, where I/T2 shows a peak. 
The value of I/T2 at the peak (Fig. |J) is also of the same 
orders of magnitude as A, consistent with our model. 

If the correlation time continuously grows at lower 
temperatures, we expect the EFG fluctuations to freeze 
eventually resulting in broadening of the NQR spectrum. 
This is indeed observed as shown in Fig. |l|. The width 
(FWHM) of 63 Cu NQR spectrum is 400 kHz at 300 K, 
while it increases rapidly below 20 K and reaches 800 kHz 
at 4.2 K. If we assume two independent contributions 
to the width, one due to T-independent inhomogeneity 
which is represented by the spectrum at 300 K and the 
other due to freezing of EFG fluctuations at low tem- 
peratures, and square of the widths from these sources 
add to make up the observed value, the FWHM due to 
EFG freezing at 4.2 K is estimated to be 570 kHz. This 
value agrees with the second moment A/2n — 170 kHz 
deduced from the 1 /Ti data, which translates to FWHM 
of 400 kHz for a Gaussian distribution. 

The results of 1/Ti, 1/T 2 , and the NQR spectra alto- 
gether provide convincing evidence for slow EFG fluctu- 
ations, which freeze randomly at low temperatures. Such 
fluctuations must be caused by motion of either electronic 
or ionic charges. Let us summarize the prominent fea- 
tures that characterize the fluctuations, irrespective of 
their origin. First, the dynamics is extremely slow, with 
1/t c ranging from 10 6 to 10 9 sec -1 in the temperature 
range 50 - 150 K corresponding to the anomalous nuclear 
relaxation. It is much too slow to be attributed to the 
transport of individual electrons along the chain. Second 
is the glassy and random nature. Gradual slowing down 
of the fluctuations without a well defined critical temper- 
ature is evidenced by the different peak temperatures of 
1/Ti and I/T2. Spatial randomness with no signature of 
a long range order breaking translational symmetry is in- 
dicated by the broadening of the NQR spectrum with no 
line splitting nor fine structure. These arc in sharp con- 
trast to what have been observed in charge-density-wave 
(CDW) materials such as NbSe3 p5| , where peak tem- 
peratures for 1/Ti and 1/T 2 coincide. Third, the ampli- 
tude of the fluctuating EFG is rather small, of the order 
of a few hundred kHz. It is known for high-T c cuprates 
that NQR frequency changes with hole concentration ap- 
proximately at a rate 20 - 30 MHz/hole f|, Hi- There- 
fore, if the NQR anomalies in Prl24 are due to electronic 
origin, only a minor fraction of the holes are relevant. 
The charge distribution at low temperatures must still 
be largely uniform and only partial random freezing is 
possible with amplitude of the order of 1 % of holes per 
site. 

Although we do not have conclusive evidence concern- 
ing the origin of the anomalous EFG fluctuations, follow- 
ing observations suggest that collective electron-lattice 
coupled motion is important. The peak in 1/Ti and 
the subsequent spectral change at lower temperatures de- 
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scribed here are very similar to what were observed in the 
lightly hole-doped two-leg ladder compound Sr24Cu2404i 
p8j. This material is an insulator at low temperatures 
and there are evidences for charge order in the ladder 
planes from both the frequency dependent conductivity 
p9[ |3(| and the development of fine structure in the NQR 
spectrum ]2§| ] . Thus the peak in 1 /T\ in this material is 
most likely caused by collective fluctuations of electronic 
charge. The similarity of the NQR results suggests that 
this is also the case for Prl24. However, Prl24 is a good 
metal and dc-transport measurements indicate no sign of 
static charge order. 

The optical conductivity o~(uj) by Takenaka et al [fl0|| , 
on the other hand, clearly shows dynamical correlation 
towards charge order instability, as mentioned earlier. 
The spectrum of a(cu) along the chain direction splits 
into the Drude like zero energy mode with only 2 % of 
the total spectral weight and a gapped finite energy mode 
centered at 40 meV. The two-component structure resem- 
bles that observed in the organic Bechgaard salts Q 
and suggests dynamic short range correlation for charge 
disproportionation. The frequency dependence of high 
energy part of <t{lo) was analyzed in the framework of 
the Tomonaga-Luttinger liquid, yielding a small value of 
the charge correlation exponent K p ~ 0.24, which implies 
strong repulsive interaction |[o| . These results altogether 
point out that Prl24 is in close proximity to charge or- 
der, therefore, charge freezing may occur near impurities 
or can be triggered by a coupling to the lattice. 

No structural transition has been reported in Prf 24 so 
far. However, recent neutron diffraction results ]3l| indi- 
cates sudden change of Ba position along the c-direction 
by about 0.05 A near f 60 K. The Raman spectra of Ba- 
phonon mode also show large frequency shift and sud- 
den narrowing near the same temperature . Although 
we notice a slight shift of the NQR frequency by about 
100 kHz (Fig. |l]) from 300 K to 100 K, which may be a 
consequence of movement of Ba, uniform change of ionic 
positions do not cause broadening of the spectrum. It is 
not likely that the static and dynamic NQR anomalies de- 
scribed here are merely lattice effects, since they are not 
observed in the isostructural Y124. On the other hand, 
when such double- well type ionic instability is coupled to 
an electron system close to charge order, we suspect that 
a static partial charge freezing can occur. 

In summary, the pronounced peaks in 1/Ti and I/T2 
at different temperatures and broadening of the NQR 
spectrum at lower temperatures for the chain Cu sites 
in PrBa2Cu408 provide evidence for slow fluctuations of 
electric field gradient. Their correlation time is extremely 
slow and show glassy freezing at low temperatures. The 
amplitude of the fluctuations and concomitant freezing 
is about two orders of magnitude smaller than what we 
expect when the whole carrier in the CU2O2 chains par- 
ticipate in charge order. 
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